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We have performed first-principle calculations of the electronic structure and exchange couplings 
for the nanotube compound Na2Vs07 using the LDA+U approach. Our results show that while 
the intra-ring exchange interactions are mainly antiferromagnetic, the inter-ring couplings are ferro- 
magnetic. We argue that this is a consequence of the strong hybridization between filled and vacant 
3d vanadium orbitals due to the low symmetry of Na2Vs07, which results into strong - and often 
dominant - ferromagnetic contributions to the total exchange interaction between vanadium atoms. 
A comparison with results of previous works is included. 

PACS numbers: 73.22.-f, 75.10.Hk 



I. INTRODUCTION 

The values of the exchange interactions between the 
magnetic moments carried by metallic atoms constitute 
a crucial piece of information for the theoretical inves- 
tigation of the magnetic properties (magnetization, sus- 
ceptibility, specific heat, etc.) of low-dimensional quan- 
tum S=i systems, for instance vanadates^ 2 * 3 - There are 
several ways to obtain this information. One of them 
is based on the estimation of the exchange integrals 
from the hopping parameters obtained by fitting the 
band-structure with an extended tight-binding model^ or 
through a downfolding 5 and projection^ procedure. The 
first-principle method for the calculation of the exchange 
interaction parameters was proposed by Lichtenstein et 
al. 7 Within this scheme the exchange interaction parame- 
ters are determined by calculating the second variation of 
the total energy S 2 E with respect to small deviations of 
the magnetic moments from the collinear magnetic con- 
figuration. 

The first theoretical investigations!^ of the magnetic 
couplings of Na2V307, a Mott insulator with the rather 
exotic topology of a nanotube, were based on the former 
method. They relied on different estimates of the hop- 
ping integrals and lead to conflicting results. Whangbo 
and Koo£ have employed a spin dimer analysis in order to 
estimate the relative magnitudes of the spin exchange in- 
teraction parameters expected for Na2V307. They have 
concluded that Na 2 V307 could be described by six helical 
spin chains and should show a gap in the spin excitations. 

Dasgupta et al& have performed the first ab-initio mi- 
croscopic analysis of the electronic and magnetic prop- 
erties of Na2Va07 using linear augmented plane wave 
(LAPW) method^ Since, due to the geometry, the 3d 
xy bands of the vanadium atoms are well separated from 
other states, these authors have derived a V 3d xy or- 
bital model Hamiltonian using the downfolding proce- 
dure within the framework of the N-th order muffin- 
tin orbital method. 11 Based on the calculated hopping 
parameters between xy orbitals, they have estimated 



the exchange interactions between vanadium atoms in 
Na2Va07 using standard superexchange theory and pro- 
posed an effective Heisenberg model. The calculated 
magnetic susceptibility as a function of temperature de- 
termined by exact diagonalization of finite clusters is in 
good agreement with experimental data. 

In this paper we report a detailed investigation of 
the isotropic exchange interactions in Na2Vs07 based 
on a first-principle LDA+U calculation of its electronic 
and magnetic structure. The exchange interactions be- 
tween vanadium atoms were calculated using the Green- 
function method developed by Lichtenstein et alH In 
order to provide a microscopic explanation of the result- 
ing exchange interactions, which are qualitatively differ- 
ent from previous estimates, we have used an extended 
Kugel-Khomskii model that includes the hybridization 
between filled and vacant orbitals using the hopping pa- 
rameters calculated by the projection procedure. We con- 
clude that estimating the exchange interactions of a sys- 
tem with a very exotic geometry such as Na2V3 07 within 
simple superexchange between xy orbitals with on-site 
Coulomb interaction is not possible. 

The paper is organized as follows. In Section II we 
shortly describe the crystal structure of Na2Va07 and 
present the results of the LDA calculation. In section III, 
we present the results of the LDA+U calculation, and we 
discuss the origin of the exchange interaction parameters. 
In Section IV, we briefly summarize our results. 



II. LDA CALCULATION 

A simplified view of the crystal structure of Na2V307 
is presented in Fig.l. The main building blocks of the 
crystal structure are distorted vanadium-oxygen pyra- 
mids. Vanadium atoms are roughly in the center of a 
pyramid which is built by five oxygen atoms. There are 
three types of vanadium atoms V1-V2-V3 which form a 
basic structure unit. Three of these units connected by 
edges are needed to form a ring of the tube. The lowest 



2 



energy orbital is the V 3d orbital of xy symmetry^ (us- 
ing a convention where the axes of the coordinate system 
are directed towards the oxygen ions in the plane) , which 
is the orbital whose lobes point away from the oxygen. 
More detailed information about the crystal structure of 
Na2Va07 can be found in Refd 




FIG. 1: (Color online) Simplified crystal structure and inter- 
action paths between vanadium atoms in Na2Vs07. Equiva- 
lent interaction paths have the same color. The tube structure 
is made by putting M and N into contact. 

The electronic structure calculation of Na2V3C"7 was 
performed using the Tight Binding Linear-Mufnn-Tin- 
Orbital Atomic Sphere Approximation (TB-LMTO- 
ASA) method in terms of the conventional local-density 
approximation^ The band structure of Na2Va07 ob- 
tained from LDA calculations is presented in Fig. 2. 
There are 18 bands near the Fermi level which have xy 
symmetry (since there are 18 vanadium atoms in the unit 
cell) and are separated from the other V 3d bands as well 
as from the 2p oxygen bands. The corresponding partial 
densities of states are shown in Fig. 3. Six bands have 
an energy significantly lower than the others, which are 
located near the Fermi level. Two of them have an en- 
ergy near -0.7 eV, while the other four are located be- 
tween -0.3 and -0.4 eV. Such a band location of active 
orbitals (occupied by electrons) is unusual in comparison 




FIG. 2: Band structure of Na2Vs07 near the Fermi level (0 
eV). 



with the band structure of other S=^ transition metal 
compounds. 15 These bands are in good agreement with 
those presented in Ref. 0- 

With one subset of energy bands well separated from 
the rest, one can expect that a tight-binding model with 
a single a;j/-orbital per V site (i.e. a 18-band model) will 
provide a good approximation to the full band structure 
close to the Fermi level. For this we express the LDA 
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FIG. 3: (Color online) Partial density of 3d vanadium states 
obtained from LDA calculations. The solid blue and the 
dashed black lines are the density of states of Sd^ and the 
total density of 3d states, respectively. 

band structure results in terms of the low-energy model 
Hamiltonian 



Htb — > t 



(1) 



where i and j are site indexes, tij is the effective hopping 
integral between sites i and j, and af a is the operator 
which creates an electron with a spin a in a state with a 
Wannier function centered at the zth atom. In order to 
obtain the hopping parameters tij , we have employed the 
projection procedure^ According to this procedure, the 
hopping integrals in the Wannier function basis [W- 71 ) 
(m is orbital index, i is site index), which is defined as 
the Fourier transform of a certain linear combination of 
Bloch functions \tpnk) (n is the band index and k is the 
wave-vector in reciprocal space) can be written in the 
general form: 



(wr\{Y.\^nm^n k \)\wf) 



(2) 

where e n (k) is the eigenvalue in the basis of linearized 
muffin-tin orbitals. The calculated inter-ring and intra- 
ring hopping parameters between 3d xy orbitals of vana- 
dium atoms are presented in Table I. These values are in 
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TABLE I: Calculated hopping integrals between xy orbitals 
of vanadium atoms for the one-orbital model and estimated 

exchange interaction parameters J?? (J*? = — jj — , U = 
2.3 eV) of Na 2 V 3 7 (in meV). 



i-j 




\i 




(Vl-1) 


-(V2-9) 


-169 


49.67 


(vi-i; 


-(V3-16) 


-145 


36.56 


(Vl-1) 


-(V2-8) 


-135 


31.70 


(V2-7) 


-(V3-18) 


-138 


33.12 


(V2-9) 


-(V3-18) 


-118 


24.22 


(V2-9) 


-(V3-15) 


-36 


2.25 


(Vl-1) 


-(V2-10) 


-32 


1.78 


(Vl-1) 


-(V3-15) 


-6 


0.06 


(Vl-2) 


-(V3-15) 


-15 


0.39 



excellent agreement with the results of the downfolding 
procedure. 9 One can see that intra-ring hopping integrals 
are in the range -120 meV to -170 meV. The inter- ring 
hopping integrals between occupied xy orbitals of vana- 
dium atoms are smaller than the intra-ring ones. 

Having set up the picture of the dynamics of individual 
electrons we shall now introduce their on-site Coulomb 
interaction U to obtain the one-band Hubbard model 

HHub = ^ i 3 a ia a j<y + 1^ n i<y n i-a , (3) 

where = af^aia is the particle-number operator. 
Based on the obtained hopping integrals, we have esti- 
mated the exchange interaction parameters using the ex- 

x 4(t lcv ) 2 

pression = — — (see Table I) , which appears upon 
mapping the one-band Hubbard model onto the Heisen- 
berg model in the limit U » t, and which only describes 
the antiferromagnetic superexchange process due to the 
hybridization between non-degenerate orbitals. The sum 
in Eq.(2) runs over all bands. Therefore in the xy-orbital 
model we take into account all interaction paths between 
xy orbitals on different vanadium atoms. 

In order to check the reliability of the mapping of the 
band structure onto the one orbital tight-binding model 
(Eq.l), we have also derived a five-orbital model de- 
scribed by a tight-binding Hamiltonian that includes five 
orbitals per V site: 

rrfull _ \ ■> ,mm'+ ( A \ 

TB ~ 2^ l iJ a im<x a im'a (4) 

m,m' 

using the projection procedure to calculate the hopping 
integrals. The results of this calculation are presented in 
Table II. One can see that the hopping parameters be- 
tween the occupied xy orbitals in the five-orbital model 
are not exactly the same as in the one-orbital model. 
This is due to the fact that the hybridization between 
filled and vacant orbitals is explicitly taken into account 
in the five-orbital model. It is also found that the sums of 
the squares of the hopping parameters between occupied 
and vacant orbitals if - are larger than those between filled 



TABLE II: Calculated hopping integrals between the 3d or- 
bitals of vanadium atoms including five orbitals per V site 
model (in eV). i™ m corresponds to the hopping process be- 
tween the filled xy orbitals of ith and jth vanadium atoms, 
while t^j™ 1 corresponds to the hopping parameters between 
filled (m) and vacant (to') orbitals. 



"j ■ i mm Ta / i mm \ 2 ~jf / < mm' \ 2 
1_ J Hj 1 ij ~ ) 1 ij ~ L^i mm I 
(m^m ) 

(Vl-l)-(V2-9) -0.166 0.028 0.055 

(V1-1)-(V3-16) -0.153 0.023 0.043 

(Vl-l)-(V2-8) -0.148 0.022 0.024 

(V2-7)-(V3-18) -0.155 0.024 0.043 

(V2-9)-(V3-18) -0.109 0.012 0.038 

(V2-9)-(V3-15) -0.042 0.002 0.036 

(V1-1)-(V2-10) -0.041 0.002 0.031 

(V1-1)-(V3-15) -0.001 0.069 

(Vl-2)-(V3-15) -0.028 0.001 0.046 



orbitals, 1°-. Therefore one can expect that the Heisen- 
berg model obtained from the one-band Hubbard model 
of Eq.(3) through perturbation theory is not sufficient in 
order to describe the magnetic couplings between vana- 
dium atoms in Na2V,307. 

III. LDA+U CALCULATION 

The main goal of the LDA+lli& method is to take 
into account the Coulomb correlations of localized 
states, which are disregarded in LDA. The former 
method provides reasonable results for strongly corre- 
lated systemsiiSiii The effective Coulomb interaction U 
and the effective intraatomic exchange J#, which repre- 
sent external parameters in a self-consistent cycle of the 
LDA+U scheme, are determined from the first-principle 
calculation by the constrained LDA method. This calcu- 
lation scheme was described elsewhere^ The Coulomb 
interaction parameter U and intraatomic exchange Jh 
have been estimated as U=2.3 eV and Ji/=0.98 eV. The 
Brillouin zone integration in our calculation has been per- 
formed in the grid generated by using (4,4,4) divisions. 

TABLE III: Results of LDA+U calculations for different types 
of magnetic configurations (see Fig. 4) . 





conf. 1 


conf. 2 


Energy Gap (eV) 


1.14 


1.18 


Magnetic moment (/is) 


0.87 


0.89 



Our LDA+U calculations show that the magnetic 
structure of Na2Va07 is frustrated with respect to mag- 
netic couplings, i.e. we cannot find a collinear spin 
configuration which corresponds to the minimum of the 
total energy of the system. We have performed the 
LDA+U calculations for two magnetic configurations 
(Fig. 4). The results of these calculations are presented 
in Table III. One can see that in both magnetic con- 
figurations Na2Vs07 is an insulator with an energy gap 
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of 1 eV. The calculated average value of the magnetic 
moment of vanadium atoms is about 0.88 he, which is 
in good agreement with experimentally observed spin-^ 
moment ^ 




FIG. 4: The different magnetic configurations for Na2VsC>7 
used in the LDA+U calculation. The arrows denote the di- 
rections of the magnetic moments, which lie along the z axis. 



orbital of site j is given by 

4"' = - r d£lm Gf L Af k ' G%% 

" J OO 

where gW = £ mm , A km G%™' (Af k ')*. In Eq.(5), the 
index n runs over all states. Therefore we take into ac- 
count all interaction paths between the ith and jth atoms. 

Our calculated values of the total exchange couplings 
are presented in Table IV. One can see that there is a 

TABLE IV: Calculated exchange interaction parameters of 
Na2Vs07 (in meV). Negative couplings denote ferromagnetic 
exchange interactions. J^ F is the contribution to the total 
exchange interaction coming from the coupling between filled 
orbitals of the ith and jth atoms, while j£ is the sum of 
the exchange interactions between filled and vacant orbitals 



of the ith and jth sites. 


i-j 






ij 


(Vl-l)-(V2-9) 


31.8 


54.2 


-22.4 


(V1-1)-(V3-16) 


13.4 


35.4 


-22.0 


(Vl-l)-(V2-8) 


22.3 


34.5 


-12.2 


(V2-7)-(V3-18) 


15.5 


35.6 


-20.1 


(V2-9)-(V3-18) 


-5.8 


16.3 


-22.1 


(V2-9)-(V3-15) 


-12.2 


~ 


-12.2 


(V1-1)-(V2-10) 


-14.0 


~ 


-14.0 


(V1-1)-(V3-15) 


-36.9 


0.4 


-37.3 


(Vl-2)-(V3-15) 


-27.7 


0.6 


-28.3 



Following Lichtenstein et ali, we determine the ex- 
change interaction parameter between vanadium atoms 
via the second variation of the total energy with re- 
spect to small deviations of the magnetic moments from 
the collinear magnetic configuration. The exchange in- 
teraction parameters Jy of the spin Hamiltonian H = 
^2i>j Jij ^iSj, where S=^, can be written in the follow- 
ing formZi2S: 

T I 7 t \ ^ / a 'mm' sym m \m"m!" /~ivn"'m\ 

J — OO , 

m" ,m"' 

where the on-site potential A™ m ' = H™ 11 ' ~ H™f and 
the Green function is calculated in the following way 

mn * (h.\ n m'n I jL\ 

Ggy'(e) = E : En ~ ^ 

k,n a 

Here c^ n is a component of the nth eigenstate, and E™ 
is the corresponding eigenvalue. In order to define the 
contributions to the total exchange interactions between 
vanadium atoms from different pairs of 3d orbitals, we 
diagonalize the on-site potential in the following way: 
A mm' = A rnk £kk (A km ')*, where A kk is the on-site po- 
tential which is diagonal in orbital space and the matrix 
A defines the basis where A is diagonal. The exchange 
interaction between the kth orbital of site i and the A/th 



qualitative and quantitative disagreement between the 
exchange interactions deduced from the hopping pa- 
rameters (Table I) and the total exchange interaction in- 
tegrals Jij obtained by the Green function method (Table 

4 (t !cy ) 2 

IV). Therefore the simple model estimation J.^ = — ^ — , 
which was used in RefH does not work in the case of 
Na2V30y. One has to take into account the low symme- 
try of this compound. We have calculated the different 
contributions to the intra-ring and inter-ring exchange 
couplings. It was found that the antiferromagnetic ex- 
change interaction between occupied orbitals jfj F (Table 
IV) is in good agreement with J^ v deduced from the hop- 
ping integrals (Table I). But there are also strong ferro- 
magnetic contributions which are due to the hybridiza- 
tion between occupied (by electrons) and empty 3d or- 
bitals of the vanadium atoms. For vanadium atoms which 
belong to different rings, one can see that the antiferro- 
magnetic contribution due to the hybridization between 
occupied xy orbitals is approximately zero (Table IV). 
Thus, the main contribution to the total exchange inter- 
action is given by hybridization between occupied and 
empty orbitals. A similar situation where the hybridiza- 
tion between filled and empty orbitals plays an important 
role was encountered for the layered vanadates CaV205, 
MgV 2 5 , CaV 3 7 , and CaV 4 0<^. 

Let us first give an illustrative geometrical explanation 
of the obtained exchange interaction picture. In order 
to do so one can consider the simple two pyramid S=^ 
model presented in Fig. 5 with the assumption that the 
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xy orbital of vanadium atom is occupied by one electron. 
The pyramid axes are collinear in the case of a perfect 




FIG. 5: Two pyramid xy model. The top and bottom figures 
correspond to perfect and distorted structures, respectively. 
The oxygen atoms are located at the edges of pyramid. 

structure, therefore only the a-overlap between the V sub- 
orbitals is non-zero. Thus the total exchange interaction 
is purely antiferromagnetic. The deviation from perfect 
geometry resulting from the rotation of a pyramid results 
in a decrease of cr-overlap between active orbitals and an 
increase of those between occupied and vacant orbitals. 
For instance, the overlap between the vacant 3z 2 — r 2 
(right pyramid) and the filled xy orbital (left pyramid) 
results in an exchange interaction of ferromagnetic na- 
ture. 

The resulting picture of magnetic interactions (Ta- 
ble IV) can be explained numerically within the Kugel- 
Khomskii modeliSi According to their investigation, for 
real materials it is necessary to consider all the 3d lev- 
els and the intra-atomic exchange interaction J# which 
leads to Hund's rule. We modify the Kugel-Khomskii 
model by introducing the splitting between the occupied 



TABLE V: Calculated crystal-field splittings between the oc- 
cupied (xy) and the vacant orbitals (yz, xz, Sz 2 -^, x 2 -y 2 ) of 
vanadium atoms obtained by using the projection procedure 
(ineV). 









A 3z 2_ r 2 




VI 


1.10 


1.71 


3.05 


2.55 


V2 


1.39 


1.08 


2.96 


2.41 


V3 


1.18 


1.60 


2.95 


2.27 



orbital (denoted as 1) and other vacant orbitals 

H = V t mm 'a + a , I V Umm ' n n . , , 

11 ~ / j L ij u ima u ']rn' a i / J „ ' Hma < W a' 



mm' 
(m^m ) 



taa 

mm 



x (1 - <Wi'<W) + / ](! — 5 lm i) (g + A m /) njm'c, 

i m' a 

^ jmm 

— ^ (1 — 3mm') 

iaa' 
mm 

X ( a tma a i"mv' a tm'a' a im'& + a tma a i™-' <? a tma' a im'a' ) ,(6) 

where £o is the energy of the filled orbital and (eo + A m >) 
is the energy of the m'th vacant orbital. For simplicity 
we assume that £o — and we neglect the fact that dif- 
ferent empty orbitals have different splittings, i.e. we 
assume that A m / = A (for vacant orbitals). It is also as- 
sumed that the Coulomb repulsion and the intraatomic 
exchange interation do not depend on the particular or- 



bital, i.e. U I; 



U„ 



U and 



3 H . The total 



exchange interaction can be expressed in terms of the sum 
of squares of the hopping integrals I?- and 1^ presented 
in Table II, the splittings between filled and empty or- 
bitals of vanadium atoms A, the on-site Coulomb U, and 
the intraatomic exchange interactions J# in the following 
form: 



J. - - J AF 



J„ 



FM 



4I a - 

i j 

U 



(U + A)(U + A-J h )' 



(7) 



The first term describes the antiferromagnetic coupling 
due to hybridization between active orbitals (occupied 
by electrons) at the zth and jth sites. The second term 
is the sum of the exchange interactions between vacant 
and filled orbitals. It is of ferromagnetic nature. 

The values of the crystal-field splittings of different 
types of vanadium atoms obtained by using the projec- 
tion procedure are presented in Table V. If we take for the 
the splitting A the average value of the crystal field split- 
tings, namely 2 eV, the ferromagnetic contributions to 



TABLE VI: Calculated contributions to the total exchange 
interactions between vanadium atoms obtained by using ex- 
tended Kugel-Khomskii model (in meV). 



i-j 




J Ah 


J ij 


(Vl-1) 


-(V2-9) 


48.7 


-26.2 


(Vl-1) 


-(V3-16) 


40 


-20.5 


(Vl-1) 


-(V2-8) 


38.2 


-11.4 


(V2-7) 


-(V3-18) 


41.7 


-20.5 


(V2-9) 


-(V3-18) 


20.8 


-18.1 


(V2-9) 


-(V3-15) 


3.4 


-17.2 


(Vl-1) 


-(V2-10) 


3.4 


-14.7 


(Vl-1) 


-(V3-15) 





-32.8 


(Vl-2) 


-(V3-15) 


1.7 


-21.9 
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the total exchange interactions estimated within the ex- 
tended Kugcl-Khomskii model are about 2 times smaller 
than those obtained using the Green function method 
(presented in Table IV) . A physically more appealing ap- 
proximation is to use in the model calculation the small- 
est value of the crystal field splitting, which corresponds 
to the first excitation energy for an electron between oc- 
cupied and vacant orbitals. In our case this splitting is 
about 1.1 eV. The calculated contributions to the total 
exchange interactions between vanadium atoms for A = 
1.1 eV are presented in Table VI. One can see that we 
can qualitatively reproduce the anti- and ferromagnetic 
contributions to the total exchange interactions between 
vanadium atoms using this simple model except for bond 
(V2-9)-(V3-18). In order to obtain full qualitative agree- 
ment with the Green function method results one should 
use A = 0.8 eV, which is slightly less than the crystal- 
field splitting calculated by the projection procedure. 

IV. CONCLUSION 

We have investigated the electronic structure and 
the microscopic origin of the exchange interactions of 
Na2V307. Based on our LDA+U results we predict 
that the value of the energy gap in electronic excitation 
spectra in this system is about 1 eV. The different con- 
tributions to the exchange interaction parameters have 
been calculated using the first-principle Green-function 
method proposed by A.I. Lichtenstein. On the basis of 



the obtained results we argue that it is not correct to 
assume that the exchange interaction which is due to hy- 
bridization between active orbitals is a good approxima- 
tion to the total exchange interaction between vanadium 
atoms for Na2Vs07. 

According to the present results, the main interactions 
are antiferromagnetic along the ring, ferromagnetic be- 
tween the rings, and of the same order of magnitude. 
The resulting Heisenberg model is strongly frustrated, 
with no obvious starting point for perturbation theory, 
hence no simple way to predict its properties. Such an 
analysis is left for future investigation. 
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